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Abstract 

Using an analytical model for the Cosmic Microwave Background anisotropics 
produced by textures, we compute the resulting collapsed three-point correla- 
tion function and the rms expected value due to the cosmic variance. We apply 
our calculations to the COBE-DMH experiment and test the consistency of 
the model with the observational results. We also show that an experiment 
with smaller angular resolution can put bounds to the model parameters. 
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I. INTRODUCTION 



An important problem in cosmology is that of finding out the source of the density pertur- 
bations leading to the large scale structure formation. Competing models are inflationary 
scenarios and topological defects, these latter formed as the consequence of a symmetry 
breaking phase transition in the early universe. 

A possible way of trying to discriminate between models is given by the anisotropics 
in the Cosmic Microwave Background (CMB) radiation predicted by them. In inflationary 
models, large scale anisotropics arise from the Sachs- Wolfe effect [0 and are proportional 
to the fluctuations of the gravitational potential on the last-scattering surface. These are 
generated as a result of quantum fluctuations of scalar fields during inflation. The resulting 
anisotropics are consistent with the analysis of the two-point correlation function of the two- 
year COBE-DMR data p[, which fixes the amplitude of the power spectrum and constrains 
the spectral index, provided that the couplings of the inflaton field are very weak. As a 
consequence, the CMB anisotropics produced in inflationary models are nearly Gaussian 
distributed. 

Defect-induced large scale structure formation has been well studied recently [^|-|^ and its 
predictions confront successfully against a large bulk of observational data. The amplitude of 
the anisotropics, which can be determined using the COBE^DMR data, is directly related to 
the symmetry breaking energy scale and corresponds to a Grand Unification Theory (GUT) 
phase transition. However, the production of anisotropy maps from numerical simulations 
of the field equations turns out to be a fabulous task, and therefore some analytical insight is 
desirable. This is the reason why people turned to consider analytical (although simplified) 
models for CMB anisotropics generated by defects, wherein the physics is more transparent 
and computations can be pursued straightforwardly [^j-^]. Although it is clear that full 
range numerical simulations will have the last word regarding these and other observable 
signatures of defects, in the meantime progress with these simplified models can hint some 
of their main features. 

CMB maps contain more information than the one that can be extracted by studying 
just the two-point correlation function. Recently, Hinshaw et al. analysed the three- 
point function of the two-year COBE-DMR data and found a non-vanishing signal, although 
consistent with the level of cosmic variance associated to a Gaussian process. The predictions 
for the three-point CMB function in inflationary models have recently been worked out in 
detail. The contribution coming from the non-linearities in the inflaton evolution and in 
the relation between the inflaton fluctuations and the CMB resulting anisotropics have been 



computed in Refs. [|ri|-|T3[. Also the post recombination mildly non-linear evolution of the 



perturbations gives rise to a non-vanishing contribution through the Rees-Sciama effect 
||l^HT6|. Both contributions result to be more than three orders of magnitude smaller than 
the associated cosmic variance, and are thus consistent with the COBE^DMR data. It is 
interesting to know whether the predictions of topological defect models are also consistent 
with the COBE^DMR three-point function data or not. As defects are the typical example 
of non-Gaussian perturbation sources, the resulting three-point function and its cosmic 
variance can easily differ from the Gaussian ones. In this paper we apply a recently proposed 
analytical model of anisotropics produced by textures to study the predictions for the 
three-point correlation function and cosmic variance of the CMB and we compare this 
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theoretical band with the COBE-DMR data. 

As an example that this statistics can be a useful test for models of structure formation, 
let us note that, based on the analysis of the three-point function of the COBE^DMR data, 
Hinshaw et al. were able to rule out an inflationary model for primordial isocurvature 
baryon (PIB) fluctuations proposed by Yamamoto and Sasaki |T^. In this model the pre- 
dicted mean skewness vanishes and the rms skewness (cosmic variance) results several orders 
of magnitude below the standard adiabatic models one, and the authors suggest that the 
rms value of the full three-point function would be tiny as well. This would make the theo- 
retical error bars for the three-point function predicted by this PIB model smaller that the 
actual data from the maps, hence running into conflict with the result of the COBE-DMR 
analysis. 

Likewise, it is interesting to check whether the non-Gaussian signatures predicted by the 
analytical texture model are in agreement with the analyses of the maps, and if not what 
the constraints are. This is the main aim of the present paper. 



II. THE COLLAPSED THREE POINT CORRELATION FUNCTION 

The three-point correlation function for points at three arbitrary angular separations a, 
(5 and 7 is given by the average product of temperature fluctuations in all possible three 
directions with those angular separations among them |jl2|. In this paper, for simplicity. 



we will restrict ourselves to the collapsed case, corresponding to the choice a = [3 and 



7 = 0, that is one of the cases analysed for the COBE-DMR data |ll8|JlO| (the other is the 
equilateral one, a = /3 = 7). The collapsed three-point correlation function of the CMB is 
given by 

C,{a) ^J^J ^AT(7i)AT2(72)5(7i ■ 72 - cos a). (1) 

For a = 0, we recover the well-known expression for the skewness, 6*3(0). By expanding the 
temperature fluctuations in spherical harmonics ^(x, 7) = Y.e,m'^Ti^)^r^i'l)^ can write 
the collapsed three-point function as 

where Tq = 2.726±0.01K is the mean temperature of the CMB radiation [IS], Wi represents 
the window function of the particular experiment and we follow the notation in Ref. []T5|, 
deflning 

Ki7r ^ J dn^y^mrmr^'f) ^ (3) 

which have a simple expression in terms of Clebsh-Gordan coefficients. 

Predictions from different models usually come as expressions for the ensemble average 
of the angular bispectrum (o-^^a^^^a^^), from which we obtain the predicted {C^{a)). This 
corresponds to the mean value expected in an ensemble of realizations. However, as we 
can observe just one particular realization, we have to take into account the spread of the 
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distribution of the three-point function values when comparing a model prediction with 
the observational results. This is the well-known cosmic variance problem We can 

estimate the range of expected values about the mean by the rms dispersion 



cv 



» = {Cl{a)) - {C,{a)y 



(4) 



III. CMB ANISOTROPIES FROM TEXTURES 

CMB anisotropics in defect theories are produced during the photon travel all the way 
from the last scattering surface to here. The computation is quite involved because it 
requires simulations following the defects, matter and photon evolution over a large number 
of expansion times @-|6;22]. In the case of textures, the non-linear evolution of the scalar 



field responsible for texture formation is the source of the CMB anisotropics. As long as 
the characteristic length scale of these defects is larger than the Hubble radius the 
gradients of the field are frozen in and cannot affect the surrounding matter components. 
However, when the size of the texture becomes comparable to (this latter grows linearly 
with cosmic time), microphysical processes can 'push' the global fields so that they acquire 
trivial field configurations, reduce their energy (which is radiated away as Goldstone bosons) 
and, in those cases where the field was wounded in knots, they unwind. In this process 
perturbations in the spacetime metric are generated, and hence also the photon geodesies 
are affected. These unwinding events are mainly characterised by one length scale, that 
of the Hubble radius at the moment of collapse, implying that the resulting perturbations 
arise at a fixed rate per Hubble volume and Hubble time. This scaling property means that 
the defect network looks statistically the same at any time, once its characteristic length is 
normalised to H~^. Then, by dividing the sky cone into cells corresponding to expansion 
times and Hubble volumes, one can restrict to the study of the CMB pattern induced in 
Hubble-size boxes during one Hubble time. Such studies have been performed by Borrill 
et al. 1^ and Durrer et al. |2^. In the latter paper, a 'scaling-spot-throwing' process 



was implemented numerically, with CMB spots derived from a self-similar and spherically 
symmetric model of texture collapse [p5| , ^ . 

An analytical approach to compute the multipole coefficients the two-point function, 
and the cosmic variances for CMB anisotropics arising in texture models have recently been 
proposed by Magueijo |Q. The model allows to study texture-induced spots of arbitrary 
shapes and relies on a hand-full of 'free' parameters to be specified by numerical simulations, 
like the number density of spots, z/, the scaling size, dg-, and the brightness factor of the 
particular spot, a^, which tells us about its temperature relative to the mean sky temperature 
(in particular, its sign indicates whether it is a hot or a cold spot). 

Let us start by discussing the spots statistics. Texture configurations giving rise to spots 
in the CMB are assumed to arise with a constant probability per Hubble volume and Hubble 
time. By dividing the whole volume in boxes of surface dS and thickness d^ we may write 
the volume probability density of spot generation in the sky as follows 

dV = -^dSdidt (5) 
H 
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FIGURES 




Observer 
(today) 



FIG. 1. Texture-induced hot and cold spots produced on different Hubble volumes and Hubble 
times, from the last scattering surface to today (see text). Dashed arrows indicate three possible 
directions on the microwave sky. The angle a corresponds to that in the collapsed three-point 
correlation function C-iia)- 

where v is the mean number of CMB spots expected to be produced in a Hubble volume 
and in a Hubble time. Now, we are interested in a two dimensional distribution of spots 
(since observations map the two dimensional sky sphere), and thus wc may integrate the 
thickness dZ in an interval of order to get the surface probability density 

dP ~ -^dSdt . (6) 

Taking into account the expansion of the universe we may express dS = r'^{t)dfl, with 
r(t) = 3t((tQ/ty^^ — 1), where t denotes proper time and to stands for its value today. It is 
convenient to change the time variable to y{t) = log2{to/t), which measures how many times 
the Hubble radius has doubled since time t up to now (see Fig. 1). For a redshift zis ~ 1400 
at last scattering we have yis ~ log2[(1400)^/^] ~ 16. In terms of y the surface probabihty 
density may be cast as 

dP = N{y)dydn, (7) 

with 

iVfe) = -2ilM?l(2»/'-l)\ (8) 

Within the model, CMB anisotropics arise from texture-induced spots in the sky. Thus, 
we may express the anisotropics as the superposition of the contribution coming from all 
the individual spots Sk produced from yig up to now. 
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AT 

— = Y,akSk{9k,y). (9) 

k 

In this expression, describes the brightness of the hot/cold k-th spot and is interpreted 
as a random variable, which characteristic values have to be extracted from numerical sim- 



ulations as those of Ref. ||23|. Sk{Ok,y) is the characteristic shape of the spots produced at 
time y, where 6k is the angle in the sky measured with respect to the center of the spot. 

A spot appearing at time y has typically a size O'^^y) ^ ds 9^°^{y), with 9^°^{y) the 
angular size of the horizon at y, and where 

• (dsH-^\ . ( 0.54 \ 

e {y) = arcsm j = arcsm [i^^j^ j • (10) 

As spot anisotropies are generated by causal seeds, their angular size cannot exceed the 
size of the horizon at that time and thus dg < 1. The scaling hypothesis implies that the 
profiles satisfy Sk{9k,y) = S{9k/9%y)). 

Combining the expression for AT/T in Eq. (|^) and its expansion in spherical harmonics, 
we easily find the expression for the multipole coefficients a™ in terms of the brightness and 
the profile of the spots 

aT = j:<^'^Si{y)Yr{%) , (11) 

k 

where 

Siiy) = 27r J'^ Sk{9k,y)P,icos9k)dcos9k . (12) 
We can now compute the angular spectrum predicted within this analytical model 

^ 7^^^ E {\aT\') = j^rf^^ E j:{akak>SiSl)Yr{lk)Yr{lk') , (13) 

+ ^) m=-e l-^^ + ^) m=-£ kk' 

where (■) denotes ensemble average. Any two different texture-spots are uncorrelated, hence 
we have 

{akak,SiSi,) = 6kk'al{si)' . (14) 

Thus, 

Ci=^T.-i{siy . (15) 

k 

The sum over all the spots contributing to the AT/T is performed by integrating this 
expression with the probability function in Eq. (^, we obtain 

Ce = To^(a^)Xf , (16) 

where we defined 

I'--'-^ jdyN{y)S'^{y)...S'-{y). (17) 

(a^) in Eq. (0) is the mean squared value of the spot brightness, to be computed from the 
distribution {ak} resulting from simulations. 
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IV. TEXTURE THREE POINT FUNCTION 



Let us now compute the angular bispectrum. Following the same steps as in the previous 
section we may express it as 



1 £2 ^3 ay ™ 1 "12 "13 



(18) 



This time the angular integral in the probability (Eq. (|^) is the responsible for the ap- 
pearance of the factor H^J^I™^^. Replacing this expression in Eq. (j|) we obtain the mean 
collapsed three-point function predicted by this model 



(C3(«)) 



(47r)2 



^ (2£i + l)(24 + l)(2£3 + l)P£i(cosa)VV,,VV,,VV,3X3^^^^^«J-,,,,,3 , (19) 

£1/2/3 



are just the squares of the the 3 j -symbols |27 



where g^^^^s = \o y 

The other quantity of interest is the expected rms dispersion about the three-point 
function due to the cosmic variance, given by Eq. (^). Computing the ensemble average of 
the combination of six a™'s, (a™^ . . . a^®), and replacing it into the expression of (T^y(a), 
we obtain after a lengthy calculation 



CV 



[a] 



2 E P£i(«)(P£i(«) + P£2(«) + ^^3(«)) -^^3 n XfW|(2£ + 1) 

£=£1 



{a^){a')Tl 



(47r)4 

r£i£i-j-£2£3£5( 



E 

^1/2/3/5,- 
•2 



Pl{a) + AP,,{a)P^,{a) + 4P,,(a)P4(a) 



^j.i.ij.2.3.5.6 >V/^W,,W,3W,,W,, J-,,,,,3^,,,3,,(2£i + 1)(24 + 1)(24 + 1)(24 + 1)(24 + 1) 



{am 

(47r)4 



E 

^1/2/3/5/ 



{a')Tl 



+^ E 

£1/2/3/4/5/6 



(47r)5 



Pl{a) + 4P,,(a)P,,(a) + 4P,,(a)P,,(a) 

•£i£2£3 j^i44 • 

P,,(a)P,,(a) W,,W,,W,3W,,W,3W,,X^^^^^^^*^^^« 



x>V|^>V,,>V,3W,,W,,J^^'^^'^^J^^'^^'^«.F,,,,,3 J-,,,3,,(2£i + 1)(2£2 + 1)(24 + 1)(24 + 1)(24 + 1) 



X J-£,£2£3-^£4£5£6(24 + l){'^^2 + 1)(2£3 + l)(2£4 + 1)(24 + 1)(24 + 1) 



(20) 



We can estimate the range for the amplitude of the three-point correlation function 
predicted by the model by {Cz{a)) ± acv{.Oi). 

Let us note that a brightness distribution {ofc} symmetric in hot and cold spots im- 
plies (a^) = and therefore yields a vanishing three-point function. This is precisely what 
happens for spots generated by spherically symmetric self-similar (SSSS) texture unwind- 
ing events [^. However, the properties of the SSSS solution are not characteristic of spot 
arising from more realistic random configurations. Borrill et al. [^] have shown that ran- 
domly generated texture field configurations produce anisotropy patterns with very different 
properties to the exact SSSS solution. Furthermore, spots generated from concentrations of 
energy gradients which do not lead to unwinding events can still produce anisotropics very 
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similar to those generated by unwindings. The peak anisotropy of the random configura- 
tions turns out to be 20 to 40 % smaller than the SSSS solution. Moreover, their results 
also suggest an asymmetry between maxima (amax) and minima (amin) of the peaks for all 
studied configurations other than the SSSS, although with large error bars. They justify 
this result as due to the fact that, for unwinding events, the minima are generated earlier 
in the evolution (photons climbing out of the collapsing texture) than the maxima (photons 
falling in the collapsing texture), and thus the field correlations are stronger for the maxima, 
which enhance the anisotropies. 



V. NUMERICAL RESULTS 

Up to this point we have left unspecified the profile of the spots, the distribution of the 
brightness factor {a^}, the typical size of the spots dg and the number density of spots z/. We 
have to fix them in order to obtain a quantitative estimate of the three-point function and the 



cosmic variance. We will essentially rely on the results of the simulations of Ref. |2^, where 
all these issues are addressed to some extent. These simulations indicate that the size of the 
spots are approximately 10% of the Hubble radius when generated, hence we take the scaling 
size dg ~ 0.1. There is no particular form of the spot profiles obtained in the simulations. 
However, a Gaussian function, as given by the profile S{9/9s{y)) oc exp(— 6*^/6'^ (y)) looks as 
a good approximation to the bell-shaped spots obtained in Ref. and we will adopt it here 
for the numerical computations. Another interesting result obtained in Ref. I^S] regards the 



number of spots produced: whereas only 4 out of 100 simulations showed unwinding events, 
the number density of nonunwinding events (which as we mentioned above also lead to spots) 
is much greater, typically of the order of one such event per simulation. This implies a mean 
number of CMB spots per Hubble volume and Hubble time i/ ~ 1. 

We need also the distribution of the spot brightness {a^} in order to compute the mean 
values of the powers of the random variable a, (a"), appearing in Eqs. ( |l9l) and (|20|) . As 
we do not have at our disposal the full distribution of values obtained in the simulations 
and we just know the mean value of the maxima and the minima, we will make a crude 
estimate by considering that all the hot spots appear with the same ah > and all the cold 
spots with the same < 0. Then, all the (a") needed can readily be obtained in terms 
of (a^) and x = {a)/{\a\). From the two-point correlation function of the COBE-DMR 
data, we can normalize the amplitude of the anisotropies and fix (a^). We do this by fixing 
Cg = (8/iK)2 which, with the help of Eq. ([T|), leads to {a'^){Tof = 1.4 x lO^K^. The 
other parameter, x, measures the possible asymmetry between hot and cold spots and we 
leave it as a free parameter. Let us note that the amplitude of the spots is directly related 
to the scale of symmetry breaking. From Ref. p3[, (a^)^/^ ~ 0.7e, with e = Btt^G^q- Hence, 



the symmetry breaking scale results $0 — 1-6 x 10~^mp, typical of GUT theories. 

We may now apply our formalism to the COBE-DMR measurements. We consider a 
Gaussian window function, Wi ~ exp{—i{i + 1)(T^/2), with dispersion of the antenna beam 



profile a = 3.2° ||2^. The analytical formulae obtained in the previous section assumed 
full-sky coverage. The effect of partial coverage, due to the cut in the maps at Galactic 
latitudes \b\ < 20°, increases the cosmic variance by a factor ~ 1.56 [|T^], and we take it into 
account by multiplying acv by this factor in the numerical results. 
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FIG. 2. Collapsed three-point function ((73(0;)) (solid line) and the rms cosmic variance (grey 
band), in 10^ fiK^ units, arising from the texture-induced temperature fluctuations. The band 
includes the ~ 50% increment in acv due to the sample variance. The asymmetry parameter is 
X = 0.07 and the dipole contribution is subtracted. 

Figure 2 shows the collapsed three-point function {C^^a)) (solid line) and the grey band 
indicates the rms range of fluctuations expected from the cosmic variance. We have plotted 



the results for x = 0.07, which looks as a reasonable value according to Ref. ||23|. It is easy 
to find how these results change for other values of the parameter x, as the central curve 
(Csi^a)) oc X, while the variance acv is approximately constant for small values of the x 
parameter. The results do not include the dipole {imin = 2), to match the COBE analysis. 
These results are to be compared with the data in Ref. [0. For all reasonable values of the 
parameter x, the data falls well within the grey band, and thus there is good agreement with 
the observations. The cosmic variance band coming from the texture model is significantly 
larger than that coming from Gaussian distributed anisotropies It is more than 

a factor 10 larger for angles smaller than 20° and a factor 5 larger for larger angles. An 
increment in the cosmic variance has also been pointed out in Ref. for the two-point 
correlation function, and has more extensively been studied in Ref. [^. The fact that the 



range of expected values for the three-point correlation function predicted by inflation is 
included into that predicted by textures for all the angles and that the data points fall 
within them makes it impossible to draw conclusions favouring one of the models, although 
one may think that the fact that the data covers the full range of the band predicted by 
inflation means that the observed sky is a typical realization of the ensemble, while it is a 
less probable realization in the texture model as all the data fall close to the origin and the 
wide range of values allowed by the cosmic variance is nearly void. 

The largest contribution to the cosmic variance comes from the small values of i. Thus, 
the situation may improve if one subtracts the lower order multipoles contribution, as in the 



10 analysis of Ref. 
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FIG. 3. Expected range of values for the three-point correlation function for the texture model 
with asymmetry parameter x = 0.4 (grey band, as in Fig. 2). Also shown is the expected range 
for inflationary models (black band). Both bands include the ^ 50% increase in acv due to the 
sample variance. All multipoles up to ^ = 9 have been subtracted. The right panel shows a zoomed 
fraction of the same plot. 

In Figures 3 a-b, we show by a grey band the expected range of values for the three- 
point correlation function for the texture model (as in Fig. 2) and by a black band the 
expected range for inflationary models. The bands do not superpose each other for some 
ranges of values of the separation angle for the value of x = 0.4 considered, what means that 
measurements in that range can distinguish among the models. The value of the parameter 



X considered is quite large, much in excess of that suggested by Ref. |]23|, but we have 
chosen it to show an example with a noticeable effect. We have checked that the non- 
superposition of the bands appears for \x\ > 0.33. The comparison with the data is less 
clear for this the two-year COBE results have large instrumental noise (in excess of 

the Gaussian cosmic variance) for this subtraction scheme [0. However, with the four-year 
data, the instrumental noise is expected to diminish under the cosmic variance, and it would 
be possible to distinguish among the models, at least for textures with |x| > 0.33. On the 



other hand, already from the present data of Ref. [|Tn[ for small and very large angles, we 
can say that models with x < —0.4 are quite disfavoured. 

An experiment probing smaller angular scales than COBE should thus be more appro- 
priate to test non-Gaussian features in texture models. As an example we compute the 
predictions for a three-beam subtraction scheme experiment with window function at zero- 
lag We = (1.5 - 2P^(cos^) + 0.5P£(cos2^))V2exp(-£(£ + l)aV2), where a = 0.64° is the 
beam width and 6 = 2.57° is the chopping angle. This window function is peaked at £ ~ 70 
and the range of multipoles that significantly contribute to the three-point function is from 
£ = 10 to £ = 100. Hence we are still probing large enough scales and our results are not 
strongly affected by the microphysics of the last scattering surface, as it has recently been 
shown that the scale for the appearance of the so-called Doppler peaks within texture mod- 
els is shifted to larger ts than in the standard adiabatic case and their contribution turns 



10 



out not to be dominant for this range of £'s []3T| , |32| . 

We obtain for the skewness S = (C^lO)) = (4.75 ± 1.94) x lO'^/iK^, where the error band 
stands for the associated cosmic variance o"^y(0), for a value of the asymmetry parameter 
X = 0.07. For comparison, the Gaussian adiabatic prediction is S'causs = 0.06 x lO^/iK^. 
Thus, in this case even for reasonably small values of the asymmetry parameter x one such 
experiment can in principle distinguish between inflation and texture predictions, and thus 
put stronger constraints on the model parameters. 



VI. DISCUSSION 

We used an analytical model of the CMB anisotropies produced by textures to estimate 
the three-point correlation function and its cosmic variance. In spite of its simplicity, we 
expect that it can give a reasonable estimation of the results that would be obtained in 
numerical simulations. The idea of this analytical model ||^ is that the texture evolution 
generates hot and cold spots in the sky at a constant rate per Hubble volume and Hubble 
time. Due to the scaling property of textures evolution the typical size of the spots is a 
constant fraction of the angle subtended by the Hubble radius at the time they are pro- 
duced. The spot shapes are given by a profile function times a brightness factor. These free 
parameters have to be determined from the results of simulations. Here we have fixed the 
number density of spots produced per Hubble volume and Hubble time, the angular scale 
of the spots and the spots profile so that they are consistent with the results of Ref. |^ . 



The computations of the three-point function and the cosmic variance are fundamentally 
sensitive to the distribution of the spots brightness. A distribution symmetric in hot and 
cold spots gives rise to a vanishing three-point function, while a non-symmetric distribution 
yields a non-vanishing one. For definiteness, we parameterized the distribution with a very 
simple ansatz with just two free parameters, namely its dispersion, (a^), and a parameter 
measuring the asymmetry of the distribution, x. {p?) was fixed using the measured ampli- 
tude of the two-point correlation function, x instead was left free and encodes the strongest 
dependence of our results on the model parameters. As discussed above, numerical simula- 
tions for random initial conditions point out to a non-vanishing asymmetry, although with 
a small value for x (and large error bars). It is straightforward to extend the calculations to 
another brightness distribution. 

The application of this analysis to the COBE-DMR experimental setting shows a very 
large cosmic variance band, much larger than that associated to a Gaussian process. For 
all reasonable values of x, this band encloses both the Gaussian one and the COBE data, 
and thus no conclusions favouring infiationary or texture theories can be obtained. The 
predictions of both models are consistent with the three-point function of the COBE^DMR 
maps. On the other hand when all the multipoles up to £ = 9 are subtracted, the situation 
looks a bit more promising and we saw that for large values of with the four years 
COBE data we could distinguish among texture and infiationary models. However, these 
large values of x seem to be larger than those expected from simulations. Finally, we showed 
that a experiment probing smaller angular scales should be able to distinguish between both 
classes of models, even for small values of the modulus of the parameter x. This looks as an 
encouraging result in the task of testing theories of the primordial fiuctuations. 
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